MULTIPLE ORTHOGONAL POLYNOMIALS AND 
A COUNTEREXAMPLE TO GAUDIN BETHE ANSATZ 

CONJECTURE 

E. MUKHIN AND A. VARCHENKO 



Abstract. Jacobi polynomials are polynomials whose zeros form the unique 
solution of the Bethe Ansatz equation associated with two sl2 irreducible mod- 
ules. We study sequences of r polynomials whose zeros form the unique so- 
lution of the Bethe Ansatz equation associated with two highest weight si r +i 
irreducible modules, with the restriction that the highest weight of one of the 
modules is a multiple of the first fundamental weight. 

We describe the recursion which can be used to compute these polynomials. 
Moreover, we show that the first polynomial in the sequence coincides with 
the Jacobi-Pineiro multiple orthogonal polynomial and others are given by 
Wronskian type determinants of Jacobi-Pineiro polynomials. 

As a byproduct we describe a counterexample to the Bethe Ansatz Conjec- 
ture for the Gaudin model. 



1. Introduction 

1.1. Content of the paper. It is well known that zeros of the classical Jacobi 
orthogonal polynomial with parameters a, (3 satisfy a system of algebraic equations, 
which is known as the Bethe Ansatz equation of the Gaudin model associated to 
sh and two irreducible modules with highest weights —a — 1, —/3 — 1, see [Hi], |V2| . 

In this paper we study the Bethe Ansatz equation of the Gaudin model associated 
to sZ r +i and two irreducible modules, one of which has an arbtrary highest weight 
and the highest weight of another one is a multiple of the first fundamental weight, 
see (|6.2ll . In this case the dimensions of the spaces of singular vectors of a given 
weight are at most one-dimensional and therefore we may expect an appearance of 
some orthogonal polynomials. In this paper we show that this is indeed so. 

We write the weights asm= (mi, . . . , m r ) and (k, 0, ... , 0), see (|2.1[l . We show 
that then for each partition I = (Zi,?2 • • ■ ,lr)> where U are non-negative integers, 
h > Ijt if i > j, and for generic values of the parameters m, k, there exists a 
unique solution of the Bethe Ansatz equation. We consider polynomials yi, . . . , y r , 
deg yi = l{, whose zeros form such a solution. 

We call parameters m, fc, I consistent if mi, . . . , m r , k are non-negative integers 
and l|2.2|l is satisfied. According to |M VI) . for consistent values of parameters, poly- 
nomials y\, . . . , y r can be computed from a certain space of polynomials V(m, I, k) 
of dimension r + 1 which has only two singular points. Degrees of polynomials in 
V(m, I, k) and the exponents at the singular points are given explicitly in terms of 
the weights m, k and the partition I, see Lemma \2. 21 
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We construct the spaces V(m,l,k) recursively using explicit linear differential 
operators, see Corollary 13.61 and discover that V(m,l,k) has a basis written in 
terms of multiple orthogonal polynomials called the Jacobi-Piheiro polynomials, 
see Lemma T4.4I and Proposition 14.61 

The Jacobi-Piheiro polynomial [F] is the unique monic polynomial of degree l\ 
whose coefficients are rational functions of m, k and which is orthogonal to functions 

1 „, J.1— fo — 1 ~— m2— 1 ~— TU2 „~rn 2 +h-l3-2 

-L ^ *Xj a ■ ■ B j *£j j *£j y *Aj a ■ ■ B - *£j B « ■ • a 

s / S / 

^^^^^^^^^ 

ii— £2 h — h 

x -T,i^2 m i- r + 1 . m 4 -r+i r 

s — ' 

with respect to the scalar product given by 

(f(x),g(x))= f 1 f(x)g(x)(x-l)- k - 1 x- m ^ 1 dx, 
Jo 

if nii, k are all negative real numbers. 

It follows that the polynomials yi , . . . , y r are Wronskians of the Jacobi-Piheiro 
polynomials up to an explicit factor. In particular y\ is exactly the Jacobi-Piheiro 
polynomial. Using properties of multiple orthogonal polynomials we are able to 
prove that for general values of parameters m, k the corresponding Bethe vector is 
non-zero, see Theorem 16. 81 

Multiple orthogonal polynomials can be computed rather explicitly and there- 
fore we have explicit formulas for y\. Similar computation allows us to obtain an 
explicit formula for y r . Using these formulas, we find that for special values of 
parameters, the sl% Bethe Ansatz equation has no solutions and therefore we ob- 
tain a conterexample to the standard form of the Bethe Ansatz Conjecture, see 
Proposition 16.101 

We also describe the spaces of polynomials U(m, I, k) which are dual to the 
spaces V(m, I, k). These spaces U(m, I, k) were considered in jSc] (with x replaced 
by —a;). In |Sc| the problem of finding an explicit basis was reduced to a linear 
system of equations of the size (X)i=i m i) ~ h- 

Thus, summarizing, we establish a relation between the Bethe Ansatz method 
and the theory of multiple orthogonal polynomials. We also establish recurrence 
relations between the vector spaces of polynomials of one variable, which have given 
dimension r + 1 and exactly two singular points at x = and x = 1. We start with 
a space of polynomials consisting of monomials and construct all other spaces using 
suitable first order linear polynomial differential operators. Moreover, we show that 
the differential operators form an extension of the sl r+ \ affine Weyl group. Those 
differential operators provide the Rodrigues-type explicit formulas for the bases of 
those spaces of polynomials. 

A study similar to the one conducted in this paper can be done in the trigonomet- 
ric case, related to multiple Laguerre polynomials, and also in the cases associated 
with the XXX model and the XX Z model, related to multiple little q-Jacobi 
polynomials studied in jPV| . See in |MV3) related discussions of the XXX model. 

Our paper is constructed as follows. In the rest of the introduction we state our 
main results in more detail for the cases of and sl^ and present our counterex- 
ample to the standard form of the Bethe Ansatz Conjecture. In Section|21we define 
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the spaces of polynomials V(m,l,k) and U(m,l,k). In Section |21 we construct 
those spaces recursively. In Section0]we construct the basis of V(m, I, k) in terms 
of Jacobi-Piheiro polynomials. In Section [3] we give an explicit basis of U(m,l, k). 
In Section [H] we use our results to study the Bethe Ansatz Conjecture. In Section 
0we present remarks on some integral identity and on an affine Weyl group action 
in our situation. 

We thank Ken McLaughlin for useful discussions. EM thanks Max-Planck- 
Institute fur Mathematik where part of this work was done for hospitality and 
a stimulating enviroment. 

1.2. Jacobi polynomials with negative integer parameters. Let m, I, k be 

non-negative integers such that 

(1.1) k>l>0, m>l>0. 

Let 

d 2 d 
T)(l,k) = x(x-l)—r - (kx + m(x -1)) — + l(k + m + l-l) . 
dx z dx 

Denote by V(l,k) the two-dimensional complex vector space of solutions of the 
differential equation D(l, k)f = 0. The space V(l, k) has the properties: 

(i) The space V(l, k) consists of polynomials of degrees I and k + m + 1 — I, 
l<k + m+l — l. 

(ii) There exists a polynomial in V(l, k) which has a root at of multiplicity 
m + 1. 

(iii) There exists a polynomial in V(l, k) which has a root at 1 of multiplicity 
k + 1. 

These three properties define the two-dimensional space V(l, k) uniquely. 

Example. Polynomials 1 and x m+1 form a basis in 1/(0,0). □ 

Let y(l, k; x) be the monic polynomial in V(l, k) of degree /. Then 

y(l,k;x) = P[ a ' 0) (x) , 

where a = —k — 1, j3 = —m — 1, and Pf a, ^(x) is the monic Jacobi polynomial on 
the interval [0, 1]. 

A basis in V(l, k) is formed by P/ _fe_1 '~ m_1) (x) and x m+1 Pilf _1 ' m+1) (x). 
According to the Rodrigues formula the polynomial y(l, k; x) is equal to 

dx l 

up to multiplication by a non-zero constant. 
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Consider the first order linear differential operators 



D (l,k) 


= x(x — 


^ dx 


- l(x-l) - k - 1 , 




D x {l,k) 


= x(x — 


x d 

^Tx 


- (k + m+1- l)(x - 1) 


- fc-1 


D% (!,*>) 


d 

dx 


- (k- 


fm+l-I), 




D?(l,k) 


d 

dx 


- i . 







The operators commute: 

D x (Z,fc + l) D {l,k) 
D?(l,k-1) D%(l,k) 
The operators define isomorphisms 



D (Z + l,fc+l) Dx(l,k) , 
D%(1 -l,k-l) D^(l,k) . 



D (l,k) 


V(l,k) - 


■* V(l,k 


+ 1), 


Di(l,k) 


V(l,k) - 


■* V(l + 


l,fc + l) 


D%(l,k) 


V(l,k) - 


■* V(l, k 


-1), 


D\{l,k) 


V(l,k) - 


+ V(l- 


l,fe-l) 



if the preimage and image spaces are defined, that is if the corresponding parameters 
satisfy (|f .ff) . In particular, the composition 

D x (l - 1, k - 1) • • • Dx(0, k - I) £> (0, k — l — 1) ■■■ D (0, 0) 

defines an isomorphism V^(0, 0) — > 1^(7, fc). 
We have 

, , , , fc + m + 2 - 21 

(1.2) y(l,k+l;x) = 



y(l + l,k + l;x) 



(k + m, + 2-l)(k + l-l) 

D-y{l,k) y(l,k;x) 



D (l,k) y(l,k;x) 



k + m+1 -21 



For any i, denote by W(f\, ■ . . , /») the Wronskian of functions fx, ■ ■ ■ ,f% of x. 
The formula 



(f,9)l,k = 



W(f,g) 



f,geV(l,k) , 



x k (x - l) m ' 

defines on V(Z, fc) a non-degenerate skew-symmetric bilinear form. 

The operators D (l, fc) and Di(l, fc) are adjoint to the operators —Dq(1, fc + I) 
and —Di(l + 1, fc + I) respectively in the following sense: 

(D (l,k)f, g)i, k +i=-(f, D%(l,k+i)g)t, k , f G V(l,k),g G V(l,k+ I), 

<£>x(Z, fc)/, S>J+i,*+i - -(/. £>i + 1, fc + l)fl)i,fc, / e V(Z, fc), G V(l + 1, fc + I). 



Let t\,...,ti be the roots of the polynomial y(l, fc; x). The roots are pairwise 
distinct, different from and 1, and form a solution to the system of equations 

m fc x - 2 
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If parameters m, I, k satisfy Iji.ljl . then system H1.3(l has a unique solution t%, . . . ,ti 
up to permutations of these numbers. 

Formulas i|1.2|l give a recurrent way to generate solutions to <|1.3|) ■ 

Coefficients of the polynomial y(l, k; x) are rational functions of k, m. 
Assume that Re k < and Re m < -1. For polynomials /, g, set 

f 1 f( x ) g( x ) 
(/- 9)k = J q xm+2 {1 _ x)k+ i dx . 

The operators Dq(1, k) and Di(l, k) are adjoint to the operators —Dq(1, k + 1) and 
—Di(l + 1, k + 1) respectively in the following sense: 

(D (l,k)f, g) k+1 = -(/, D^(l,k + l)g) k , 
{Dx(l, k)f, g) k+ i = -(/, D^(l + l,k + l)g) k , 

We have 

(y(l,k;x), x n ) k = , 1 < n< I , 

and 

i m + 1 - i 

{y(l,k;x), y(l,k;x)) k = l\ ]J — — — — ^—^ (1,1)*-; , 

where 

r(-m- 1) T(-k + l) 



(1, l) fc _i 



r(-fc-m- 1 + 



1.3. Three-dimensional spaces of polynomials with two singular points. 

The parameters m — (mi, 7712), £ = (ii, £2) €E Z> , fc S Z>o are called consistent if 

(1.4) k > h > l 2 > , mi > h-h, m 2 > h . 

Let m — (mi,m 2 ), £ = (h,l 2 ), k be consistent. Let 

d 3 d 2 d 

D(l,k) = x(x-l)—r - ((2 mi +m 2 )(x -1) + 2kx)—r + v x - h v , 

dx^ dx* dx 

where 

v-i = -l\ + l 2 {m 2 + 1) +h(k + mi + 1 - h + h) + 

, . , . mi (m-r + m? + 1) fc(fc + 1) 
(m!+A:)(mi+m2 + A: + l) — £-J_^' v 



x x — 1 

(fc + + l 2 (m 2 + 1) + h(m x +k + l-k+ l 2 )) _ 

x{x — 1) 

h(h ~ h +m 1 + k + 1)(2 -l 2 + m 1 + k + m 2 ) 



Denote by V(l,k) the three-dimensional complex vector space of solutions of the 
differential equation D(l, k)f — 0. The space V(l, k) has the properties: 

(i) The space V(l, k) consists of polynomials of degrees l\ , mi + 1 + k + 1 2 — l\ , 
and mi + m 2 + 2 + k — l 2l where 

l\ < mi + l + fc + i2 — ^1 < mi + m-2 + 2 + k — l 2 . 
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(ii) There exist three polynomials in V(l, k) which have a root at of multi- 
plicities 0, mi + 1, mi + m 2 + 2, respectively. 

(iii) There exist three polynomials in V(l, k) which have a root at 1 of multi- 
plicities 0, k + 1, k + 2, respectively. 

These three properties define the three-dimensional space V(l,k) uniquely. 

Example. Polynomials 1, x mi+1 , and x mi+m2+2 form a basis in V(0, 0), where 
0=(0,0). □ 

Set 

e (l, k) = li, ei(i, k) = mi + 1 + k + l 2 — h, e 2 (Z, k) = mi + m 2 + 2 + k - l 2 . 
Set 

T^x) = x mi (x~l) k , T 2 (x) = x m * . 

Let /o,/i be the monic polynomials in V(l,k) of degrees e (l,k), ei(l,k), respec- 
tively. Set 

(1.5) Vi(!,k;x) = f (x) , y 2 {l,k;x) = rV^TTT^vi r ( \ ' 

(eo(t, k) - ei(t, k)) ii(x) 

These are monic polynomials of degrees l\ and h, respectively. 
Let 

d 3 d 2 d 

T> v (l,k) = x(x-l)— - ((2m 2 + m 1 )(x-i) + kx)— + u x — + u , 

where 

Mi = h(h -h + m 2 + l) + (h + m 2 + l)(m 2 +m 1 +k + 2-l 2 )- 

m 2 (m 1 + m 2 + 1) 



M 



—2 — 2m 2 — mi — k — 
h(h -l 2 + m 2 + 1)(— 2 - m 2 - mi - fc + /i) 



a; 

Denote by ?7(Z, fc) the three-dimensional complex vector space of solutions of the 
differential equation D V (Z, fc)/ = 0. The space U(l, k) has the properties: 

(i) The space U(l, k) consists of polynomials of degrees l 2 , m 2 + 1 + l\ — l 2 , 
and mi + m 2 + 2 + k — Zi, where 

Z 2 < m 2 + 1 + Zi - Z 2 < mi + m 2 + 2 + fc - Zi . 

(ii) There exist three polynomials in U(l, k) which have a root at of multi- 
plicities 0, m 2 + 1, mi + m 2 + 2, respectively. 

(iii) There exist three polynomials in U(l, k) which have a root at 1 of multi- 
plicities 0, 1, k + 2, respectively. 

These three properties define the three-dimensional space U(l,k) uniquely. 

Example. Polynomials 1, x m2+1 , and x mi+m2+2 form a basis in {7(0, 0). □ 
Set 

eQ (i, /c) = Z 2 , ei(Z, k) = m 2 + 1 + Zi - l 2 , e 2 (I, k) = mi + m 2 + 2 + k — h. 
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For fx, f 2 ,f 3 eV(l,k), set 

w U f n _ W^/1,/2) w t rf . n W{h,h) 

For given fx,f 2 jfsj the function Wy-(/i, /2) is a polynomial, and the function 
Wl{fi,f 2 ,h) is constant. The map (fi,f 2 ,fo) ^ ^(/i,/ 2 ,/ 3 ) defines a vol- 
ume form on V(l, k). 

For 31,32,53 G */(*,*), set 

Wir (31 > 52) = 7f , ^£7(31,32,33) = ^2^ • 

J 2 J- 2 ±X 

For given 31,32,33, the function W^(<7i, 32) is a polynomial, and the function 

^4(31,32,33) is constant. The map (31,32,33) >-» ^£7(31,32,33) defines a volume 
form on U(l,k). 
We have 

[/(Z,fc) = {W^(/i,/ a ) I /i,/a G V(/,fc)} , 

V(l,k) = {^(31,32) I 31,32 G 17(1,*)} . 

In particular, the polynomial y 2 (l,k;x) defined in 1)1.5(1 is the unique monic poly- 
nomial in U(l, k) of degree l 2 . 

The spaces V(l, k) and (7(Z, k) are dual with respect to the following pairing: 

(/, 9}i,k = Wl(f, fx, h) , if g = Wl{hJ 2 ) . 
In particular, we have 

( yi{l,k;x), y 2 (l,k;x) }j ifc = . 

Consider the first order linear differential operators 

Di(l,k) = x(x-l)^- - ei(l,k)(x-l) - fc-1, 

A V M) = x± - 4_&,k) , 

ax 

i = 0, 1, 2. There are linear relations: 

A l2 {l,k) D (l,k) + A 20 (l,k) Di(i,fc) + D 2 (i,k) = 0, 

4 2 (i,*) £»o(^*) + A w (l,k) DX{l,k) + A 01 (l,k) D y 2 (l,k) = 0, 
where Aij(l, k) = ei(l, k) — ej(l, k). Set 

1 = (0,0) , li = (1,0) , 1 2 = (1,1) . 
The operators commute: 

Di(l + lj,k+ 1) Dj(l,k) = Dj(l + li,k + l) Di(l,k) , 
Dy(l-lj,k-l) D^(l,k) = Dj(l-li,k-l) D^(l,k) 

for all i, j. 

The operators define isomorphisms 

D t (l,k) : V(l,k) -» + , 
£>/(*,*) : 17(1, A) — »■ 17(1 — li, A — 1) , 
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for i = 0,1,2, if the preimage and image spaces are defined, that is if the corre- 
sponding parameters satisfy (|1.4(l . In particular, the composition 

D 2 ((Ii-l,/ 2 -l),A-l) ••• D 2 ((h-l 2 ,0),k-l 2 ) Di^Ii-Ia-l.Oj.fc-Ia-l) 
••• Di((0,0) s fe-ii)I>o((0,0),fc-li-l) ••• D ((0,0),0) 

defines an isomorphism V(0, 0) — > V(l, k), and the composition 

£ o v ((0,0),l) ••■ D^((0 7 0),k-h) D^((l,0),k-h + l) 

••• D^((h-l 2 ,0),k-l 2 ) D%((h -l 2 + l,l),k-l 2 + l) ■■■ D y 2 ((hM),k) 

defines an isomorphism U(l,k) — > 17(0,0). 
Set Iq = k, I3 — 0, and 



Mhk) = n 



k + J2)=i TUj — + 1 



o k + J2]=i nij — h+ i + l — li+l 



i+l 



j4i (i, k) = 2l\ — l 2 — k — m\ — 1, A 2 (l, k) = l± + l 2 — k — mi — m 2 — 2, 
A^(l,k) = h + l 2 -k-mi-m 2 -2, Ai(l,k) = 2l 2 -h-m 2 -l, 



Ej= 3 -i m 3 - h + i 



Sj=3-i TO j _ h + i + 1 - h-i + h-i 
Then for i = 0,1, 2, we have 

^(Z + lj.fc + lja;) = Di(l,k) yi(l,k;x) , 
Ay{l,k)y 2 {l-li,k-l;x) = D^(l,k)y 2 (l,k;x) . 

Here is a Rodrigues type formula for y±(l, k;x): 



rai+m2+2 



(X-1) 



k+l 



dx J 



„— ni2+l-2 — 1 



ix-Ia 



h— 1 2 — mi— lf x _ fe— 1 

= c yi(l,k;x) , 



where c is a non-zero constant. If 

h h 

yi(l,k;x) = ^2 a n (l,k)(x -l) n , y 2 (l,k;x) = b n {l,k)x n 

n=0 n=0 

with a h (l,k) = 1, a_i(Z, fc) = a h+1 (l,k) = 0, and b h (l,k) = 1, &; 3+1 (Z,fc) = 0, 
then 



a„(Z + li,fc+l) 

(1.6) 6 n (Z + l i5 A; + l) 
In particular, 



fc + l + n n-l-e^Z,*;) 

+ 7-77— tt a n _i(i,«) , 



4(1,*) 
A^(l + l u k + l) 



Ai(l,k) 
b n (l, k) . 



b n {l,k) = 
where Iq = k. 



h-n-l 2 

n n 



T,l=3-jm s -l 2 + i + j 



i=o r- 



1 E^3-, TO s - h + i +j + 1 + h-j - h-j 
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We have 

A 12 (l,k) A (l,k) yi (l + l ,k+l;x) + 

A 20 (l,k) Ax(l,k) yi(Z + li,fe + l;x) + 

A i(i,fe) A 2 (l,k) yi(l + l 2 ,k + l;x) = 

A 12 (l, fc) Alii, fc) y 2 (l -l Q ,k— 1; x) + 

A w (l,k) A\(l,k) y 2 (l~l 1 ,k-l;x) + 

A 01 (l,k) A%{l,k) y 2 (l-l 2 ,k-l;x) = 0, 

and 

yi(I,M) # 0, yi(l,k;l) + 0, |ft,(I,fc;l) ± 0. 
Example. If Z = (2, 1), then 

yi ((2,l),fc;x) = (x-i) 2 + 

(fc- l)(2fc + 2mi +m 2 ) , . fc(fc - 1) 

-(a;-l) + 



(mj + fc — 1)(toi + 77i2 + k) (m 2 + mi + k)(m± + fc — 1) ' 

W ((2,d,m) = « - ? ^T 2 ^ 

(m 2 + 2)(fc + mi + m 2 ) 
In particular, if (2m i + m 2 ) 2 + k(Amx — m 2 ) = 0. Then 
yi((2,l),*;a:) = lfc((2, 1), k; x) 2 . 

□ 

For i = 0, 1, 2, the operator -Dj(i, fc) is adjoint to the operator — D^(l + l i: k+ 1) 
in the following sense: 

(Di(l,k)f, g)i +li , k+ i = - (/, A v (^ + l,^ + l)5)j,fc , 
if / S V(l, fc), and geU(l + l h k + 1). 

Let ti, . . . , i/j be the roots of the polynomial fc; x). Let si, . . . , s; 2 be the 
roots of the polynomial y 2 (l,k;x). Assume that yi(l,k;x) and y 2 (l,k;x) do not 
have multiple roots, then the roots satisfy the system of equations 

(1.7) 

mi fc x - 2 \ - 1 

— + + H = ' i = l,--.,h, 

ti ti 1 .... ti tj ti Sj 



3i J^i 

m 2 x - , x • . 

^ ^_ ' 1 = 1, ■••,<2 

3, 3i^i 3 



2 ^ 1 



If parameters m,l,k are consistent, then system (|1.3fl has at most one solution 
ti, . . . , t; 15 Sx, . . . , si 2 up to permutations of the first and second groups of these 
numbers. 

If parameters I are fixed, then for almost all m, k such that m, Z, fc are consistent, 
the polynomials y\(l, k;x) and y 2 (l, k\x) do not have multiple roots. 
Formulas l|1.6|l give a recurrent way to generate solutions to l|1.7|) . 
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Coefficients of the polynomials yi(l, k; x) and 1/2(1, k; x) are rational functions of 
k, m. 

Assume that Re k < and Re mi + m 2 < —2. For polynomials /, g, set 

ft x _ f 1 m 

{L 9)k ~ J Q (1 _ x )k+l aX ■ 

For i = 0, 1, 2, the operator Di(l, k) is adjoint to the operator —D^(l + li, k + 1) 
in the following sense: 

(Di(l,k)f, g) k+ i = - (/, D^(l + li,k + l)g) k . 

We have 

yi(l,k;x) x n 



1 yi(l,k]x) x n 



dx = , n = 1, . . . , li , 

dx = , n = 1, . . . , h — 1% 



l x m i+ 2 (1 - x) k+1 
We also have 

(yx(l,k;x), y 2 (l,k;x)) k = C(l,k) , 

where 

r(-k + h) r(-TO! - m 2 - 2) 



= Z 2 ! 



r(— fe — mi — ui2 — 2 + Zi ) 



n 1 (m 2 - i) (mi + m 2 + 1 - i) 1 yr 1 m 2 + 2 + i 

n (m 2 + h-l 2 + l-i){k + m 1 +m 2 + 2-l 1 -i) 2 k + mi + l-h-i' 

In the next formulas we include m in the notation for yi(l,k;x) and V(Z, fc). 
Then for consistent m, I, k, the space V(m, I, k) has a basis formed by polynomials 

yi((mi,m 2 ), {h,l 2 ),k;x) , 

x mi+1 yi((-mi - 2, mi + m 2 + 1), (k + l 2 - l u l 2 ), k; x) , 

x m 1+ m 2 +2 yi (( m2j _ TO1 _ m2 _ 3)j ( fc + ; 2 _ Zl) _ fe; x ) . 

1.4. Application to the Bethe ansatz in the Gaudin model. We let g = 

n + © f) © n_ be a simple Lie algebra, a\, . . . ,a r £ f)* simple roots. 

For a g-module L and fi € ()* denote by L[/j] the weight subspace of L of weight 
fi and by SingL[/j,] the subspace of singular vectors of weight fi. 

Let n be a positive integer and A = (Ai, . . . , A„), Aj G f)*, a set of integral 
dominant weights. For a weight fi € f)*, let L M be the irreducible g-module with 
highest weight ^. Denote by La the tensor product Lai © ■ ■ • © La„ ■ 

Let z = (#1, . . . , z n ) be a point in C™ with distinct coordinates. Introduce linear 
operators Ki(z), . . . , K n (z) on La by the formula 

K i(z) = 2^ — — , i = l,...,n. 

3; 1=^1- 

The operators are called the Gaudin Hamiltonians of the Gaudin model associated 
with L\. The Hamiltonians commute. 

The problem is to diagonalize simultaneously the Hamiltonians, see B , BF , FFR 
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The Hamiltonians commute with the g- action on La- It is enough to diagonalize 
the Hamiltonians on the subspaces of singular vectors, Sing L\ [//] C La- 

The eigenvectors of the Gaudin Hamiltonians are constructed by the Bethe 
ansatz method. 

Fix /i = Yl s =l ^ s — S[=i where li, . . . , l r are non-negative integers. Set 

+ _ ,(2) (2) (r) ,(r)s 

b — yb 1 , . . . , t ;i , L 1 , . . . , L l2 , . . . , L 1 , . . . , j . 

One defines a suitable rational function w(t,z) with values in £a[m] as in |RV| . 
cf. |MV2I IR.SV| . The function is symmetric with respect to the group Sj = 

(i) 

x • • • x S/ r of permutations of coordinates t - with the same upper index. One 
considers the system of equations 

l 



(1.8) _v^^ + E Ett^ + E 4^1 = °' 

s=l <*j z s Sj s ^ik=l L j L k s, s=jtj L j L s 

where i — 1, . . . , r and j = 1, . . . , Z,. The system is symmetric with respect to Sj. 
One shows that if t° is a solution to i|1.8[l . then u>(i ,z) belongs to Sing La [fj] and 
2;) is an eigenvector of the Hamiltonians Ki(z), . . . , K n (z), see |RV| . 

This method of finding eigenvectors is called the Bethe ansatz method. System 
(|1.8|) is called the Bethe ansatz equation, the vector w(t° , z) is called a Bethe vector. 

The Bethe Ansatz Conjecture says that if Ai, . . . , A„ are integral dominant, the 
weight fj, is integral dominant, and zi,...,z n are generic, then the Bethe vectors 
form a basis in Sing La [/x]. In particular, the conjecture implies that the number 
of Sj-orbits of solutions to 1)1. 8|) is not less than the dimension of Sing La [m]. 

Assume that g = s?3, n = 2, zi = 0, z% = \. Assume that for some non- negative 
integers fc, mi, m2, the weights Ai and A2 are such that 

(Ai,ai)=mi, (Ai,a 2 ) = m 2 , (A 2 ,ai) = fc, (A 2 ,a 2 ) = 0. 

Then system l)1.8fl becomes system (ll.7|) . 

Let = Ai + A 2 — Ziai — liOii be integral dominant, then the space of high- 
est weight vectors SingLAi ® L\ 2 [fi] is one-dimensional, if parameters m,l,k are 
consistent, and is zero-dimensional otherwise. 

It was shown in MV1 that if /i is integral dominant and the numbers t\, . . . , U x 
and Si, . . . , si 2 form a solution to l|1.7fl . then m, I, k are consistent and 

h h 

Y[(x-ti) = yi(l,k;x) , JJ(a;-Si) = y 2 {l,k;x) . 

i=l i=l 

Then we obtain, see the Example in Section ll. 31 

A counterexample to the Bethe Ansatz Conjecture. 

Let I = (2, 1). Let m,l,k be consistent and 

(1.9) (2mi+m 2 ) 2 + k{Am x ~m\) = 0, 
then 

yi ((2,l),k;x) = y 2 ((2,l),k;x) 2 , 

and roots of these two polynomials do not form a solution to p.7| ) . Hence system 
i 1. 7| ) does not have solutions in this case. Hence there is no Bethe vector to generate 
the one- dimensional space Sing (La 1 <E> La 2 )[m]- 
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For example parameters m = (2,3), k = 49, I = (2,1) are consistent and 
satisfy ljl.9jl . 

2. Spaces of polynomials 

In this section we define our main objects - the spaces of polynomials with two 
singular points 0, 1 and special exponents at 1. 

2.1. An sl r+ i lemma. Fix a natural number r and choose non- negative integers 
mi, . . . , m r - Choose a non-negative integer k. 

Consider the complex Lie algebra sl r +i- Let f) be its Cartan subalgebra, let 
a>i,...,a r E t)* be its simple roots, ( , ) the scalar product on f)* such that {on, on) = 
2. For a weight A G f)* denote by L\ the irreducible sl r +i module with highest 
weight A. 

Let the weights Ai,A2 be such that the scalar products with simple roots are 
given by 

(2.1) (Ai,Q! s ) = m s , (A 2 ,a s ) = S al k, 

where s = 1, . . . , r. Thus L\ 2 is the fc-th symmetric power of the vector represen- 
tation and Lai can be any irreducible finite-dimensional sZ r +i module. 
We describe the decomposition of the tensor product <g> L\ 2 . 

Lemma 2.1. We have an isomorphism of sl r+ \ modules: 

where the sum is over non-negative integers l\, . . . , l r such that 

(2.2) k > h > h > ■ ■ ■ > lr > 0, l s - l s +i <m s (s = 1, . . . , r). 

In particular the multiplicity of every module in the right hand side is one. □ 

We use the notation / = (li, . . . , l r ) and m = (mi, . . . , m r ). We call I a partition 
if U are integers and ^i > I2 > • ■ ■ > l r > 0. We set 

^0 = k, l r+ i = 0. 

We call parameters m, I, k consistent if all m i; ij and k are non-negative integers 
satisfying inequalities (|2.2|) . 

2.2. Definition of spaces of polynomials V(m, I, k) and U(m, I, k). Fix a posi- 
tive number r. Let V be an (r+l)-dimensional complex vector space of polynomials 
of one variable x. Assume that V has no base points, i.e. for any zeC there exists 
/ € V such that f{z) ^ 0. 

The Wronskian polynomial W(V) ofV is the Wronskian of any basis fi, . . . , f r +i 
of V: 

W(V) = W(f 1} f r+l ) = det (£r . 

The Wronskian polynomial W(V) is defined up to multiplication by a non-zero 
complex number. 

A point z E C is called a singular point of V, if z is a root of W(V). 
For any z E C there exist unique non-negative integers n\{z), . . . , n r {z) such that 
for any / £ V the order of zero of / at z belongs to the set 

r 

(2.3) { 0, 1 + ni(z), 2 + ni(z) + n 2 {z), . . . , r + ^ n »( z ) }• 

3=1 
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A point z G C is a singular point of V if and only if at least one of the numbers 
ni(z), . . . ,n r (z) is positive. The set ()2.3[l is called the set of exponents ofV at z. 
For zeCwe define a full flag in V 

5{z) = {Q = F Q (z) C Fx{z) G F 2 {z) C • • • C F r+1 {z) = V}, 
where the i-dimensional subspace F^(z) consists of all polynomials with root x = z 
of order at least r — i + 1 + X)s=i +1 n s{z). We call the flag $(z) the flag of V at z. 

Similarly there exists unique non- negative integers Ui(oo), . . . , n r (oo), I such that 
for any / G V the degree of / belongs to the set 

r 

(2.4) {I, 1 + 1 + ni(oo), 1 + 2 + m(oo) + n 2 (oo), . . . , I + r + n s (oo) } . 

3=1 

The set ()2.4|) is called the set of exponents of V at infinity. 
We define a full flag in V 

y(oo) = {0 = F (oo) G F x (oo) G F 2 (oo) c ■ • • C F r+1 (oo) = V}, 
where the i-dimensional subspace Fj(oo) consists of all polynomials of degree no 
greater than I + i — 1 + 53s=i n « (°°) ■ We call the flag ^(oo) t/ie flag of V at infinity. 
We have 

W(V) = Y[(X - (.r+l-s)n a 

and therefore the following relation 

r r 

(2.5) ^ ( r + 1 _ s ) n ^ = ( r + !) ^ + 2 ( r + 1 ~ s ) n »(°°) ■ 
zeCs=i s=i 

Given a full flag J = {0 = F c Fl c • • • C F r+i = V^} we call a basis 
{/x, . . . , /r+i} °f V compatible with flag 5F if {/i, . . . , / s } is a basis of F s for s = 
1 r 

We use the notation n(oo) = (ni(oo), . . . ,n r (oo)), n(z) = (ni(z), . . . ,n r (z)). 

Fix arbitrary sequences of non-negative integers n(z),n(oo) so that they satisfy 
relation l|2.5|) . Let A(z), A(oo) be sZ r +i integral dominant weights defined by the 
conditions 

(A(z),aj) = rii(z), (A(oo),Q!i) = n,-(oo) (« = 1, .. . ,r). 

Then it is known, see EH , MV1 , that the number of spaces V with the given 
n(z), ti(oo), I is finite and bounded from above by the multiplicity of the module 
AiY(oo) in the tensor product ® z L^ z y Moreover the number of such spaces is non- 
zero if and only if the above multiplicity is non-zero. 

In this paper we study spaces of polynomials which have only two singular points 
Z\ = and z 2 = 1. Moreover, we consider spaces only with special exponents at 1. 
Namely, we restrict our attention to two cases 

ra(l) = (Jfe, 0, . . . , 0) and n v (l) = (0, . . . , 0, k). 
For our purposes we parameterize the exponents with the above restriction using 
the parameters of Section 12.11 Given parameters m, I, k we define two sets of 
parameters n(0), n(l), I, n(oo) and n v (0), n v (l), l v , n v (co) by the formulas: 

(2.6) Uj(l) = S u k, ni(0)=mi, 

I — l\, Hj(co) = li—i + Zj-i-i — 11^ + rtii + Suk (i = 1, . . . , r) 
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and 
(2.7) 

n^(l) = 5 ir k, 7i 4 v (0) = m r+ i-i, 

l v = l r , 1^(00) = i r _, + l r +2-i — 2Z r+ i_i + m r +i-i + S ri k (i = l,...,r). 

Both n(0),n(l),/,n(oo) and n v (0), n v (l), Z v , n v (oo) satisfy relation d^75f > . 

In this parameterization the sets of exponents at infinity given by (|2.4|l take the 
form {e (m, I, k), . . . , e r (m, I, k)} and {eo(m, i, fc), . . . , e^(m, i, fc)}, where 

(2.8) ej(m, J, fe) = fc + 2^m s - + + i, 

s=l 
r 

(2.9) e^(m,l,k)= ^ m s - Z r - i+i + lr-i + h 

s— r+1 — i 

where i = 0, . . . , r. Note that for < i < j < r we have 

(2.10) ej(m,l,k) - ei(m,l,k) = e^^ra, - e^_j(m,l,k). 

Specializing the general statements to our case we obtain: 

Lemma 2.2. If parameters m,l,k are consistent then there exists a unique space 
of polynomials V(m, I, k) of dimension r + 1 im'i/i exactly two singular points and 
1 and exponents given by §2.8} and \2.4\ where n(0), n(l), n(oo) are given by §2.6} . 
If parameters m,l,k are not consistent then there exists no such a space. □ 

Lemma 2.3. If parameters m,l,k are consistent then there exists a unique space 
of polynomials U(m,l,k) of dimension r + 1 with exactly two singular points 
and 1 and exponents given by §2.8} and §2.4$ where n(0), n(l), n(oo) are equal to 
n v (0), n v (l), n v (oo) given by §2.61 . If parameters m,l 7 k are not consistent then 
there exists no such a space. □ 

Our parameters m are often fixed. We often drop the dependence on m from 
our notation and simply write V(l, k), U(l, k). 

2.3. Connection between spaces V(m, I, k) and U(m, I, k). Let m, Z, k be con- 
sistent parameters, V = V(m,l,k) and U = U(m,l,k) the corresponding spaces 
of polynomials. 

We define a tuple of polynomials T = (T\, . . . , T r ) by 

(2.11) T 1 = (x-l) k x m \ Ti=x mi (i = 2,...,r). 
Lemma 2.4. We have 

W{V) = T{T r 2 - x . . . T r , W(U) = T r r T r T Z\ ...T ± . 

□ 

For fx , . . . , fi £ V and gx, ■ ■ ■ ,gi £ U define the divided Wronskians 

Wl(fx, . . . , fi) = W(fx, . . . , fi) T}-*!*-* . . . Tr\ , 

Wl{gx, ...,<&•) = W(gx,. ..,<&) T^T^l ■ ■ ■ T r ~-i+2- 

We also define divided Wronskian Wy(Vx) (resp. Wjj(t/i)) of any subspace Vx C V 
(resp. Ux C U) as the divided Wronskian of any basis in Vx (resp. Ux). The 
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divided Wronskians W^(Vi), W}j(Ui) arc defined up to multiplication by a non- 
zero constant. 

The divided Wronskians are polynomials. 

Lemma 2.5. We have 

V = {W^{gi,...,g r ), gi,...,g r e U} , 

U = {Wl(fu---Jr), /l,.-.,/r eU} . 

Proof. The vector spaces V and {W}j{jg\, ... ,g r ), gi, ■ ■ ■ ,g r G U} have the same 
dimensions, the same singular points, and exponents, see |MV1| for details. Hence 
they coincide by Lemma \'2. 21 □ 

2.4. Examples of spaces. Let mi,...,m r be non-negative numbers. Let = 
(0, . . . , 0). Then the parameters m, 0, k — are consistent. 

Consider the space Vo = V(m, 0, 0). Let e 2 ; = ei(m, 0, 0). Then we have 

i 

e j = i + rrij , 

3=1 

where i = 0, . . . , r, see (|2.8f) . 
Lemma 2.6. We ftawe 

V" = span{l = x e °, x ei , x er }. 

□ 

Proof. The exponents of the space Vo at and infinity coincide. Thus the polyno- 
mial in Vo with the highest order of zero at has to be x £r . Similarly, a polynomial 
in Vq with the order of zero at equal to e r _i has to be a linear combination of 
x^- 1 and x £r . Continuing this argument we obtain the lemma. □ 

Fix a non-negative integer I, and let I =(/,... , I). Assume that m r > I. Then 
the parameters m : l° , k = I are consistent. 

Consider the space Ui = U(m,l°,l). Let e 4 v = e^(m,P,l). Then we have 

i 

eo = l , &i = ^(m,l°,l) = i + 2jm r+ i_j , 

3=1 

where i = 1, . . . , r, see (|2.9|) . 

Lemma 2.7. There exists a unique monic polynomial Uq(1) of degree I such that 

V V 

Ui = span{u (7), x £l , . . . , x £r }. 
Moreover, uq(1) = x l — Bix 1 ^ 1 + .... where 

s =i e o e « 1 

Proof. The last r exponents of J7; at are equal to the last r exponents of Ui at 
infinity. Thus the functions x £l , . . . , x £r belong to Ui . Therefore there exists a 
unique monic polynomial uq(1) of degree I such that the functions Uq(1), x £l , ■ ■ ■ , x £r 
form a basis of {/;. Moreover, we have 
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where c is a constant. The first non-trivial coefficient Bi of the polynomial uo(l) is 
easily computed from this formula. □ 

Corollary 2.8. We have 

C/ =span{l = x e o( m ' '°), /( mA0 », z<(™.°>°)}. 

□ 

3. Recursion for spaces of polynomials 

In this section we explain how the spaces of polynomials with different values of 
parameters are related to each other. 

3.1. Recursion of the first type. Fix a sequence m of non-negative integers. 
Let I, k be such that m, I, k are consistent parameters. 
Introduce linear first order differential operators 

d 1 

(3.1) Di(l, k)=x{x-l)—-(k + J2m s - h + h +1 + i)(x - 1) - k - 1, 

s—l 

where i = 0,1, ... ,r. Note that the coefficients of (x — 1) are exactly the degrees 
of polynomials occurring in the space V(l, k). Namely, the coefficient of (x — 1) in 
Di(l, k) is equal to ei(l, k), see l]2.8fl. 

The operators Di(l,k) are linearly dependent. Let 

(3.2) A ab = A ab (l, k) = e a {l, k) - e b (l, k). 

Lemma 3.1. The vector space spanned by operators Di(l,k), i = 0,...,r, has 
dimension two. Moreover, for 0<i<j<s<r, we have 

A js Di(l,k) + A si Dj{l,k) + A l3 D s {l,k) = 0. 

□ 

Set li = (1, . . . , 1, 0, . . . , 0) where we have i ones and r + 1 — i zeroes. 
The operators Di possess the following holonomy property. 

Lemma 3.2. For all i,j E {0, 1, . . . , r} we have 

Dj(l + l h k+l) Di(l,k) = Di(l + lj,k + l) D 3 (l,k). 

Proof. It is an easy explicit check. □ 

We give a conceptual reason for Lemma f3. 21 in Section ET31 

Theorem 3.3. Suppose m, I, k and m,l + li, k + l are consistent parameters. Then 
operator Di(l, k) maps V(l, k) to V(l + 1,, k + 1) and 

Di(l,k): V(l,k) -^V(l + li,k+ 1) 

is an isomorphism of vector spaces. 

Proof. We have 

Ker Dj(l,k) = c(x - l^x^^i^-h+h+i+J-^ 
where j = 1, . . . , r and c is an arbitrary constant. 
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Note that the kernel of Dq contains no non-trivial polynomial. Recall that 
numbers J"? 1 m s + i — 1 and 5Z* =1 m s + * are consecutive orders of zeros at x = 
which appear in V(l, k). Under our assumptions we have 



i — 1 i i 

m s + i - 1 < m s — Zi + ij+i + i — 1 < ^ ; 

s— 1 s— 1 s—1 



Therefore the intersection of Ker Di(l, k) and V(l,k) is trivial, because the orders 
of zero at x = of functions in these two spaces do not match. 

It follows that V := Di(l, k)(V(l, k)) is a space of polynomials of dimension r+1. 
Let us consider the exponents of V. 

Let fi, . . . , fr+i be a basis in V(l, k) compatible with the flag at infinity. We 
have deg f s = e s , where e s are given by l|2.8|l . Then we obviously have 

deg D l f s = deg f s + 1 (s^i), 
degA/i < deg/,. 

Also if g € V and g(l) = then the order of vanishing at x = 1 of Di(l, k)g(0) is 
at least k + 2. 

It follows from relation l|2.5|) that the above inequality is an equality and the 
exponents of V coincide with the exponents of V(l + li,k + 1). Therefore those 
two spaces coincide by the uniqueness part of Lemma 12.21 □ 

Corollary 3.4. Let m,l,k be consistent and m,l + lj,fc + 1 are not consistent. 
Then U — h+i = and 

( x -i) k+1 x^ i ^ m '- li+k + 1+i - 1 £ V{l,k). 

In particular dimZ)j(i, k)V(l, k) = r. 

Proof. If the statement of the lemma were wrong, then by the same argument as in 
the proof of Theorem ^. HI the space Di(l, k)V(l, k) would be a space of polynomials 
V(l + li,k + 1), which does not exist by Lemma \'2. 21 □ 

Lemma 3.5. Let m^l, k be consistent parameters. 

• There exist a sequence of indices ■ ■ ■ , i(k) £ {0, . . . , r} such that I = 
J^ =1 lj( s ) and m, ^s=i ■'•»(«) > 3 are consistent parameters for j = 0, . . . , k. 

• This sequence is unique up to permutation of its elements. 

• occurs k — h times and i occurs U — U+\ times [i = 1, . . . , r). 

Proof. Note that the addition of 1.; adds one to U — and leaves all other differ- 
ences lj — lj+i, j i intact. In particular the addition of lo does not change either 
of these differences. The lemma follows. □ 

Corollary 3.6. Let l,k be such that m,l,k are consistent. There exist a unique 
linear differential operator Di k of order k and of the form 

fc-i 

(3.3) D Lk := D m £l j( ,,,fe-l) ... D m (l j(1) , 1) D j(1) (0, 0) 

s=l 

with j (s) G {0,...,r}, T, k s=i 1 j(s) =1, such thatV(l,k) = A,fc V(0,0). 

Proof. The existence follows immediately from Theorem 13.31 and Lemma 13.51 The 
uniqueness follows from Lemma 13. 51 and Corollary 13. 41 □ 
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Sometimes we stress the dependence of the operator Di y k on m and write D m .i,k- 

We give some remarks about a basis in V(l, k). 

Lemma 3.7. The space V(l, k) has a basis {vq,Vi, ...,v r } such that v is a polyno- 
mial of degree l\ and Uj, i = 1, . . . , r, have the form 

mi— 

Vi = (x - l) k+1 X^=l a *3 x3 > 

3=0 

where ay are non-zero constants. 
Proof. The lemma is true for Vo, 

TCli 

3=0 

In addition if D = x(x — l)d/dx + A(x — 1) — (k + 1), then 

(3.4) D : {x-l) k+1 x a i — > {k + l + a + A) (x-l) k+2 x a . 

These two statements imply the lemma. □ 

Note that using (|3.4|l . one can write the constants ay explicitly in a factored 
form. 

Note also that Corollary 13 . 41 follows from Lemma f3. 71 

3.2. Recursion of the second type. The construction of this section is parallel 
to that of Section |3~T1 

Fix a sequence m of non-negative integers. Let l,k be such that m,l,k are 
consistent parameters. Introduce linear first order differential operators 

j r 

(3.5) DY(l,k)=x—-( Y m s -l i+1 +k + r-i), 

dx ^-^ 

s=i+l 

where i = 0, . . . , r. Note that the constant coefficients in operators DY(l,k) are 
exactly the degrees of polynomials occurring in the space U(l,k). Namely, the 
constant coefficient in DY(l,k) is exactly e^._ i {l 1 k), see l|2.9[) . 

The operators DY are linearly dependent. Recall constants A a {, given by (|3.2ll . 

Lemma 3.8. The vector space spanned by operators D^(l,k), i — 0, . . . , r, has 
dimension two. Moreover, for 0<i<j<s<r, we have 

A Ja DY(l,k) + A sl D](l,k) + AijD^(l, k) = 0. 

□ 

The operators DY possess the following holonomy property. 
Lemma 3.9. For all i, j 6 {0, 1, . . . , r} we have 

Dj(l-li,k-l) DY(l,k) = £>y(J-lj,fc-l) Dj(l,k). 

□ 

We give a conceptual reason for Lemma 1331 in Section 1X^1 
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Theorem 3.10. Suppose m,l — lj,fc — 1 and m^l^k are consistent parameters. 
Then operator DY (I, k) maps U(l, k) to U(l — lj, k — 1) and 

I)/ (I,*): U(l,k) -» U(l - U,k- 1) 

is an isomorphism of vector spaces. 

Proof. The proof is parallel to the proof of Theorem 13. 31 □ 

Corollary 3.11. Let l,k be such that parameters m,l,k are consistent. There 
exist a unique linear differential operator k of order k and of the form 

k 

Dtu := D) {i) (l m , 1) D] {2) + l j{2) ,2) ... D^ k) (£ l j(s) , k) 

s=l 

with j(s) £ {0, . . . , r}, J2 k s =i lj(s) = I, such that U(0, 0) = D^ k U(l, k). □ 

3.3. Canonical maps. Fix a sequence m of non-zero integers. Let h, fa and I2, fe 
be such that the parameters m,li,ki and m^li^k?, are consistent. 

Recall that given any space of polynomials V and zeCwe have the full flag of 
V at z, ?(V, z) = {0c Fi(V, z) C F 2 (V, z) C • • • C F r+ i(V, z) = V}, defined by the 
order of vanishing of polynomials in V at x = z. We also have the full flag of V at 
infinity, ?(V, 00) = {0 c F X (V, 00) c F 2 (V, 00) c • • • C F r+1 (F, 00) = 7}, defined 
by the degrees of polynomials in V. 

We call a linear map a : V(Zi, fci) — > ^(^2; fe) canonical if a is an isomorphism 
and it is compatible with flags at singular points and infinity in the following sense: 

a(Fi(V(l u h), 0)) - Fi(V(l 2 , A*), 0) (i = 1, . . . , r), 

a(F(V(Zi, fei), 00)) = F(F(/ 2 , k 2 ), 00) (1 = 1,..., r), 
a(F r (^(/i, fa), 1)) = F r (^(Z 2 , fa), 1). 

We call a linear map b : U(l\, ki) — > f/(Z2 ; fa) canonical if 6 is an isomorphism 
and it is compatible with flags at singular points and infinity in the following sense: 

b(Fi(U(l u fa), 0)) = Fi(U(l 2 , fa), 0) (i = 1, . . . ,r), 

b{Fi{U{l u fa), 00)) - F(t/(Z 2) fa), 00) (* = 1, . . . ,r), 
b{Fx{U{h, fa), 1)) = F X {U{1 2 , fa), 1). 

Clearly inverse of a canonical map is canonical and a composition of canonical 
maps is canonical. A canonical map multiplied by a non-zero number is canonical. 

Theorem 3.12. Let m,li, ki and m,Z 2 , fa be consistent parameters. Then there 
exist unique (up to multiplication by a non-zero number) canonical maps V(l\, fa) — > 
V(I 2 ,*2) andU{l ll k 1 )^U{l 2l k 2 ). 

Proof. We give a proof for the case of spaces V(l, k). The spaces U(l, k) are treated 
similarly. It is obviously sufficient to consider the case of (Zi,fa) = (0,0), that is 
the case of V(l\, fa) = Vq. 

We have a map Vb — > T^(Z 2 ,fa) given by the operator -D; 2 ,fc 2 , see Corollary 13. 61 
This map is clearly canonical, see the proof of Theorem 13.31 

Since we have the existence, to prove that the canonical map is unique it is 
sufficient to consider canonical maps Vb — > Vq. 
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Let a: Vq — ► Vq be a canonical map. Let 3^(0), 5(oo) be the corresponding 
flags at 0, 1 and oo. By degree reasons we have 

(3.6) Wl o (Fi(oo)) = 1 (i = l,...,r) 

(up to non-zero constants). 

For i = 1, . . . ,r + 1 consider Gi = Fi(oo) n F r +a— i(0). By the definition of 
a canonical map, a(Gi) = Gj. But because of (|3.6|) we have dimGi = 1. Let 
Gi = cfi, fi G Vq. The basis {/i, . . . , f r +i} is compatible with 5"(oo) and the basis 
{/r+ii • • • i /i} is compatible with ^(0). In the basis {/i, . . . , f r +i} the map a is 
diagonal. Let a(f s ) = fi s f s . 

Similarly for i = 1, . . . , r + 1 we have i^(oo) H = c$i for some gi G Vo- 

The basis {gi, . . . , g r +i} is compatible with 3"(oo) and the basis ■ ■ ■ ,<?i} is 

compatible with The span of 32, ■ ■ ■ , 3s is the intersection of F r (l) and F s (oo) 

and therefore it is preserved by the map a. 

Without loss of generality we can put fx = g\ = 1. It follows from Lemma \l. 61 
that without loss of generality we can assume 

It follows that if we write fi = Yl , v S ig s then va ^ and vy t ^ 0. 

Since the map a preserves the span of gi we get /ii = (12 — M- Since the map a 
preserves the span of <?2, 33 and this span does not contain g±, we get fi^ = /i. And 
so on. 

It follows that the the map a is scalar. □ 

Note that holonomy properties described in Lemmas 13.21 and 13.91 are natural to 
expect because of Theorem 13.121 since both compositions are the unique canonical 
maps when acting on our spaces of polynomials. 

3.4. Inclusion Property. Fix a sequence raofr non-negative integers. Let a se- 
quence I of r non- negative integers and a non-negative integer k be such that m, I, k 
are consistent parameters. Let s be a natural number. Let I = (Zi, . . . , l r , 0, . . . , 0) 
where we have s zeroes, and fh — (mi, . . . , m r , rii, . . . , n s ), where rij are non- 
negative integers. 

Lemma 3.13. We have 

V(m,l,k) c V(ffi,l,k). 

Proof. The lemma is true for the case of I = 0, k = 0, see Lemma l2~Bl The canonical 
isomorphisms a : V(m, 0, 0) — > V(m, I, k) and a : V(fh, 0,0) — > V(fh,l,k) are 
given by the same formula a = D m .ik — j k — a. The lemma follows. □ 

Canonical maps are compatible with the inclusion described in Lemma 13.131 
Namely, let I 1 , fc 1 be such that m, Z^fc 1 are also consistent parameters. We set 
T = (l{,...,ll,0,...,0) where we have s zeroes. We have canonical isomorphisms 

a : V(m,l,k) -> V(m, I 1 , k 1 ), a : V(m,T,k) V(fh,t, k 1 ), 
defined by 

a = D m l i k i o (Drni^y 1 , a = D^jx^ o (D^j ifc ) _1 . 

Corollary 3.14. The canonical map a is the restriction of the canonical map a to 
V(m,l,k) G V(rh,l,k). □ 
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Spaces U (m, I, k) possess similar properties. For example, the analog of Lemma 
13. 131 for spaces U(m, I, k) has the following form. 

Lemma 3.15. We have the surjective well-defined map of vector spaces of polyno- 
mials 

7r : U(fh,l,k) — > U(m,l,k), 

dx dx dx 

□ 

3.5. Duality. Fix a sequence m of non-negative integers. Let I, k be such that 
m, I, k are consistent parameters. 

We use Lemma f2. 51 to define a pairing ( , )i.k : V(l, k) eg) U(l, k) — » C as follows. 
For / G V(f, fc), G L7"(l, fc) we set 

{f,g)l,k = W^ (m {f,fx,...,f r ), if fl = W^ (I(fc) (/l,...,/r), 

see Section 6.1 in MV1 , where this pairing is studied in a more general setting. 

Lemma 3.16. The pairing ( , )i^ ■ V(l, k) <8>U(l, k) — > C is a well-defined bilinear, 
non-degenerate pairing. 

Proof. Lemma f3 . 1 61 follows from Lemma \'2. 51 □ 

Recall that the degrees of polynomials e, and ej occurring in V(l, k) and U(l,k) 
are given by (|2.8[) and 1)2. 9fl . 

Lemma 3.17. Let /, <? &e monic polynomials such that f G V(Z,A;), deg/ = e,, 
and gr G U(l,k), deg 5 = e^L^-. Then we have 

r 

(f,g)i,k= ]J ( e «- e 0- 

Proof. For a! < ■ ■ • < aj we have the formula 

W{x a \ . . . , x a i) = x^=^ ".-itf+D/ 2 Y[(a s - a t ). 

s>t 

We compute (denoting the terms of lower degrees by dots) : 
(f,g)l,k = 

W l{l,k)(f^ e ° + ■■■,■■■ ,x-^T..., ... ,x e * + ...)\ [] ( e ^ e *) 



(-iy]J(e s -e t )i [] (^-e t )\ = J] ( e - ~ ^ 



□ 



Proposition 3.18. The linear operators Di(l,k) : V(l,k) — » y(Z + lj,fc + 1) 
and —D^(l + lj, A; + 1) : Z7(Z + lj, A; + 1) — > E/(Z, fc) are adjoint. Namely for any 
f G V(Z, fc) and any g G C7(Z, fc) we /lave 

(A(/ ) fc)/,9)!+i ( ,fe+i = -(/,A y (' + ii.fc + i)flKfc- 
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Proof. The operator Di(l,k) defines a canonical map. It follows from definitions 
that the map adjoint to a canonical map is canonical. Therefore, by the uniqueness 
of the canonical map, the map defined by D^(l + 1,, k + 1) is a constant multiple 
of the map adjoint to the one defined by Di(l, k). 

We claim the constant is — 1. Let fo, ■ ■ . , f r be a basis of monic polynomials 
in V(l,k) compatible with the flag at infinity. We have deg/ s = e s and fi are 
monic. Let e s = e s + 1 — <5i. s be the degrees of polynomials in V(l + li,k + 1). 
Let g G U(l + lj, k + 1) be a monic polynomial of maximal degree. In the next 
calculation we assume that i > 0, the case of i = is similar. By Lemma 13.171 

(Di(l, k)f , g) = (e - e l )(x h+1 + ...,<?> = 

r 

= (eo - ei) Yi ( g s - eo) = (e - + 1) JJ(e s - e ). 

s, s^O s=l 

Here we denoted the terms of lower degrees by dots. Similarly we obtain 

T" 

(fo, DY(l + li,k+ l)g) - (e v r - e r v _ 4 - 1) [](e s - e ). 

s=l 

Therefore by ^111)) we have (Di(l,k)f ,g) = -(f ,D?(l + U,k + l)g) and the 
constant is —1. □ 

4. Spaces V(m, I, k) and Jacobi-Pineiro polynomials 

In this section we present an explicit basis of spaces V(m, I, k) in terms of mul- 
tiple orthogonal polynomials, called Jacobi-Piheiro polynomials. 

4.1. Rodrigues formula. Let k, 1%, . . . , l r be non-negative integers such that k > 
^i > • ■ ■ > ^r- Let mi, . . . , m r be generic complex numbers. 

For i = 0, . . . , r we define a differential operator which acts on functions of x: 

bJm, I, k) = s£i=i "*+* f ! ' +1 x li - li+1 -^=i m *-*. 

This operator is the composition of multiplication by a monomial, multiple differ- 
entiation, and another multiplication by a monomial. 

Lemma 4.1. The differential operators Di(m,l, k) commute for all values of pa- 
rameters: 

b i (m 1 ,l 1 ,ki)b j {rn2,l2,k 2 ) = bj(m 2 , h, k 2 )bi(mi, h, fa). 

□ 
Set 

(4.1) oj(x) := (x - l)- k - l x -^=^ mi - r . 

Recall operator D m i j, given by i|3.3fl and our conventions l Q = k and l r +\ = 0. 
Proposition 4.2. We have the equality of differential operators 

D m j, k = (x- l) k+1 b r {m, I, fe)D r _i(m, l,k) ... D (m : I, k) ' 



x- 1 



uj-Hx) 



a „lr-lr+i-rrw-l 



lr — lr+l—mr — 1 ™U — 1 — U m% — 1 — 1 



x h— h— mi— 1 



d x lo-h x - 1 
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Proof. For any complex a, b and non-negative integer s we have the following iden- 
tity of differential operators: 

x(x — I)— a(x — 1) — b — s } ( x(x — 1)- a(x — 1)— b — s + 1 

ax J \ ax 

x(x - 1)4- - a(x -l)-b]=(x- i ) fc + s +V- fe+1 4—tt( x - iy b x- a+b+s . 
y ' dx J dx s+1 

Indeed, both operators are of order s + 1 with the symbol (x(x — l)) s+1 d s+1 . More- 
over, it is easy to check that both operators kill the functions 

(x - l) b+l x a - b - s (i = 0,...,s). 

It follows that these two operators are the same. 

The proposition follows from the above identity. □ 

4.2. Weyl group action. Let k,h, ■ ■ . ,l r be non-negative integers satisfying the 
conditions k > l\ > • • • > l r . Let mi, . . . , m r be generic complex numbers. 

Let S r +i be the the Weyl group of sl r +i. Then § r +i is isomorphic to the group of 
permutations of the set {0, 1 . . . , r}. We denote w(i) the image of i € {0, 1 . . . , r} un- 
der permutation w. The group S r +i is generated by simple transpositions s%, . . . , s r . 
The simple transposition Sj permutes elements i and i — 1 and preserves all other 
elements of the set {0,1, ... ,r}. 

The group § r +i acts on si r +i weights by the formula 

St ■ A = A - (A + p, cti)a.i. 

The first formula in 12.1f> identifies sl r+ i weights with vectors m. Under this 
identification the Weyl group action becomes: 

Si • m — (mi, . . . , m,_2, m,_i + m< + 1, -m, - 2, m i+ i + mi + 1, m i+2 , m r ). 

We denote the result of the action of w G § r +i on m by w ■ m. 

Given k, define another action of the Weyl group § r +i on the sequences of 
numbers I by the formulas 

(si)fcZ = (h, . . . , h-i, h-i + U+i — U, k+i, ■ ■ ■ , l r ), 

where as always Iq = k and l r +\ — 0. We denote the result of the action of w G § r +i 
on I by (w) k l- 

If we consider the map C r — > C r+1 mapping I — > {Iq — l\, l\ — l 2 , ■ ■ ■ , l r — l r +i}, 
then the second Weyl group action descends to the standard permutation action of 
coordinates. 

The Weyl group acts on operators Di(m, I, k) as follows. 

Lemma 4.3. We have 

D t ( Sl -m,( Sl ) k l,k) = x- m '- l D sl(i) (m,l,k) x m '+\ 
Di(sj ■ m,(sj)kl,k) = D S](i ){m,l,k), 

where j — 2, . . . , r. □ 



24 



E. MUKHIN AND A. VARCHENKO 



4.3. Basis in V(m,l, k). Let m,l,k be consistent parameters. 
The following lemma follows from Theorem 13.31 and Lemma 12.61 

Lemma 4.4. The space V{m 1 l, k) has a basis of the form 

{v {m, I, k), «!(m, I, k)x mi+1 , v r {m, I, k)x^=^ m ' +r }, 

where Vi(m, I, k) is a monic polynomial of degree 

deg Vi(m, l,k) = k — k + U + i 

whose coefficients are rational functions of m with coefficients in Q. Moreover, the 
polynomial v%{m,l, k) is obtained via the action of operator D m _i_k as follows 

(4.2) Vi (m, l,k)=x- D m , Lk (xZU m 

□ 

The polynomials Vi are permuted by the Weyl group. 
Proposition 4.5. Let w G S r +i be an element of the Weyl group, then 

Vi(w ■ m,(w) k l,k) =v w (i)(m,l,k). 
Proof. For i = 2, . . . , r from Lemmas 14 . 31 and 14 . II we have the equality of operators 

and 

r) _— mi— In mi+1 

ly s 1 -m,(s 1 ) k l,k — LJ m,l,kX 

The proposition follows from formula (|4.2|l . □ 

Let P(m,l,k) be the monic polynomial defined by application of differential 
operators: 

P(m, I, k) = c(x - l) k+1 D r (m, I, fe)5 r _i(m, l,k) ... Di(m, I, k) ((x - l) 1 ^- 1 ) = 
cwHx) 4-r-r- " 



d X lr-lr+l dx lr + 1 ^ lr 

_! d h ~ l 

dx h ~ l i 



h— i3-"i2-l_ ^i— 2a— mi— 1^, _ j^Zi-fc-lj 



where c is a non-zero constant. 

Following |AfjV| . we call P(m,l,k) the Jacobi- Pineiro polynomial. Clearly, 
P{m, l,k) is a polynomial of degree l\ in x, whose coefficients are rational functions 
in m, k with coefficients in Q. 

Sometimes we drop the dependence on m from our notation and simply write 
P(l,k). 

Proposition 4.6. Let w G § r +i be an element of the Weyl group, such that w(0) = 
i . Then 

Vi(m, I, k) = P(w ■ m, (w)kl, k). 

Proof. In view of Lemma 14.51 it is sufficient to prove the proposition for w = id. 
We have equality of functions 

v (m,l,k) = L> m ,z,fc (1) = 

(a; - l) k+1 D r (m,l,k)D r ^x(m,l,k) . ..D Q (m,l,k) (— ^ j • 
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But 



1 \ fjk-h / k-h ' 

D (m, l,k)[ = -_- 7 )={k-h)l{x- l) h -> ! 



x—l) dx k ll \x—l / 
The proposition follows. □ 

Corollary 4.7. Let w E § r +i be such that w(0) = 0. Then 
P(m,l,k) = P(w ■ m,(w)kl,k) . 

□ 

We also note that the divided Wronskian of all elements of the basis is a constant, 
which gives an identity for Jacobi-Piiiciro polynomials. 

Corollary 4.8. We have 

W(P(m,l,k), P{ Wl ■m,{w 1 ) k l,k)x mi+1 ,...,P(w r ■ m,(w r ) k l,k) x^=^ m ' +r ) = 

const i^-i( r+1 ^'"' (x- l) rk , 

where Wi — Sj . . . S2S1. 

4.4. Properties of Jacobi-Pineiro polynomials. Jacobi-Piheiro polynomials 
are examples of multiple orthogonal polynomials, see [F], |IN| . |ABV| and refer- 
ences therein. In this section we prove some properties of Jacobi-Pineiro polyno- 
mials, which we need for applications. 

We fix a sequence m of arbitrary complex numbers. 

Define constants Ai(l, k) by 



A (i,k)= n 



k + J2t=i m s — h+s + l 



k + Z)i=i m s - h+s + l + l s +i - l s ' 

i 

Ai(l, k) = e (l, k) - ei(l, k) = h — }^m s — k + lj - 



where i = 1, 
We have 



A (i,k) = (k-h+i) Yi 



e s (l, k) - e (l, k) + l s - h+i + 1 



e s (l,k)~e (l,k) + l 
These numbers satisfy the following property. 
Lemma 4.9. We have 

^(Z + ljjfc + l) Aj(l,k) = Aj(l + li,k + l) Ai(l,k). 

□ 

The recursion for spaces of polynomials translates to the recursion for the Jacobi- 
Pineiro polynomials: 

Proposition 4.10. We have 

Di(l, k) P(l, k) = Ai{l, k) P(l + U,k + 1). 



26 



E. MUKHIN AND A. VARCHENKO 



Proof. First we suppose that parameters m, I, k and m, I + lj, k + 1 are consistent. 

By Theorem 13.31 and Proposition ^. 61 we obtain that Di(l, k)P(l, k) is a constant 
multiple of P(l + li,k + 1). Call these constants Ai(l, k). 

Recall that P(l,k) is monic. Comparing the highest degree terms we obtain 
Ai(l,k) — Ai(l,k) for all cases except for i = when the highest degree terms 
cancel and we need to consider terms of the next degree. 

By Lemma f3. 21 the numbers Aj(Z, k) satisfy the relation of Lemma f4. 91 We also 
have the initial conditions 

i (o,fc) =k + l. 

The numbers Aq(1, k) are uniquely determined by Ai(l,k), (i = 1, ...,r) and 
these properties. But numbers Ao(l,k) satisfy the same recursion. Therefore 
A (l,k) = A (l,k). 

The case of general parameters m,l,k follows by analytic continuation with 
respect to m,k. □ 

Lemmas 13.11 and Proposition 14. 1UI immediately imply the following 3-term rela- 
tion for Jacobi-Pineiro polynomials. 
Let I be a partition. 

Corollary 4.11. We have 

AiAj S P{1 + l h k + 1) + AjA is P(l + lj, k + 1) + A s Aij P{1 + 1 S) k + 1) = 0, 

where A t = A t {l,k), A tp = A tp {l,k). □ 

5. Spaces U(m,l,k) 
In this section we compute explicitly a basis in the space U(m,l, k). 

5.1. Basis in U(m, I, k). Let m, I, k be consistent parameters. Recall that we use 
the convention Iq = k. 

The following lemma follows from Theorem 13. 101 

Lemma 5.1. (Cf. |Sc|. ) The space U(m,l,k) has a basis of the form 

{u (m, I, k), u\(m, I, k)x mr+1 , . . . , u r {m, I, fc)x^i =1 mi+r } j 

where Ui(m,l, k) is a monic polynomial of degree 

degUi(m,l,k) = Z r _j - l r - i+ i 

whose coefficients are rational functions ofm,k with coefficients in Q. □ 

The polynomials Ui are permuted by the Weyl group. 

Lemma 5.2. Let w G § r +i be an element of the Weyl group, then 

u r -i(w ■ m, (w)kl, k) = u r - w (i) (m, I, k). 

Proof. We have the formula 

(5.1) u r ^(m,l,k) = c i W\f ,h,...,f i ,...,f r )x-^=r-i + ^-i j 

where Cj are non-zero constants, fi = Vi(m, I, k)x^">=i ms+l and Vi are as in Lemma 
14.41 Now the lemma follows from Lemma T4.5I □ 

Note that formula (|5.1|l determines Uj in terms of Jacobi-Pineiro polynomials. 
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5.2. Recursion and 3 term relation. Fix any sequence of r complex numbers 
m. Define constants A^(l,k) by 

r 

AY(l, k) = e^(l,k) - e^_i(l,k) = l r - ^ m s ~ h + h+i - r + i, 

s=i+l 

Z)i=i m s -l r + s 



A 



(i,fe) = -JJ_ 

s=0 ?-'i=l 



m s — l r + s + 1 + l r — s — I 



r — s+1 



where i = 0, . . . , i — 1. 

These numbers satisfy the following property. 

Lemma 5.3. We have 

^(i-l^fc-l) A){l,k) = A) {I - l u k - 1) A^{l,k). 

In addition 

A w ((k,...,k),k) = B k , 
where B k is defined by □ 

The recursion for spaces of polynomials U (m, Z, k) translates to the recursion for 
the first polynomial uq. 

Proposition 5.4. We have 

DV(l, k) u (l, k) = A^(l, k) u (l -li,k- 1). 

Proof. Similar to the proof of Proposition 14 . 1 01 □ 

Let I be a partition such that I — lj, I — Lj, I — l s are partitions. 

Corollary 5.5. We have 

A^A JS u (l - h, k - 1) + A)A is u {l - lj, k-l) + A^A tj u (l — l s ,k — T) = 0, 
where A? = A%(l,k), A tp = A tp (l,k). □ 

5.3. Explicit formulas. Proposition l5.4l allows us to compute the polynomials Ui 
explicitly. 
Write 

i r 

uo(m, I, k) = ^(-l) l Ci(m, I, k)x lr ~ l . 
Proposition 5.6. We have 

I 1 b\- f lr \TJTT T,t=r-s+l m t-lr + S+j 

Cil ' ' } \i ) ji l\ E r t =r- S+1 m t -l r + s + j + l + l r - s - l r - s+1 ■ 

Proof. We have co(m, I, k) = 1. By Proposition l5.4l we have the following relations: 

(l r — i — e^_j)ci(m, I, k) = Aj(m, I, k)a(m, I — lj, k — 1). 
So, the proposition is obtained by multiplying all the constants. □ 

Note that in fact we computed all Ui because of Lemma 15.21 
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6. Applications to the Bethe Ansatz 

We consider a special case of Bethe Ansatz equation. It has at most one solution. 
If this solution exists then it is given by zeroes of Jacobi-Pineiro polynomials and 
of Wronskians of Jacobi-Pineiro polynomials. We use properties of Jacobi-Pineiro 
polynomials to prove that for generic values of weights, this solution does exist and 
the corresponding Bethe vector is non-zero. However, we find that in special cases 
the Bethe equation has no solutions and we obtain a counterexample to the Bethe 
Ansatz Conjecture for the Gaudin model. 

6.1. Generalities of the Bethe Ansatz method. Let g = n + © f) © n_ be a 

simple Lie algebra, (., .) the Killing form on g, ct\, . . . , a r € f)* simple roots. 

Let (xi)i£i be an orthonormal basis in g, £1 = Y] ie j Xi®Xi £ j the Casimir 
element. 

For a g-module V and /j£f)* denote by V[pi] the weight subspace of V of weight 
ll and by SingV^/i] the subspace of singular vectors of weight /i, 

SingT/|/i] = { v S V | n+v = 0, hv — (li, h)v } . 

Let n be a positive integer and A = (Ai, . . . , A„), Aj 6 f)*, a set of weights. For 
a weight /i £ fj*, let L M be the irreducible g-module with highest weight /i. Denote 
by La the tensor product Lai ® • • • ® La„ ■ 

If X e End (La;), then we denote by X( l > £ End (La) the operator • • ■ ® id <g> 
X <g> id <g> ■ ■■ . If X = Y,k X k ® Y k e End (L A , cg> L Aj ), then we set X^ = 

E k 4°®n W) eEnd(L A ). ^ 

Let z = (z\, . . . , z n ) be a point in C" with distinct coordinates. Introduce linear 
operators Ki(z), . . . , K n (z) on La by the formula 

K i{z) = 2^ — - — , i = l,...,n. 

3: 3¥* 

The operators are called the Gaudin Hamiltonians of the Gaudin model associated 
with La- The Hamiltonians commute, \Ki(z), Kj(z)\ = for all 

The problem is to diagonalize simultaneously the Hamiltonians, see B , BF , FFR 

OIEnilE21E21E2] 

The Hamiltonians commute with the g-action on La, thus is enough to diago- 
nalize the Hamiltonians on the subspaces of singular vectors, Sing La C La- 

The eigenvectors of the Gaudin Hamiltonians are constructed by the Bethe 
Ansatz method. 

Fix Aqo = Y2s=i ^ s ~ li a i wnere h,- ■ - ,l r are non-negative integers. Set 

+ _ (Al) Al) p) A2) Ar) Arh 

\ b l 5 ■ ■ ■ 1 b l 1 1 b l 5 • • • ) L l 2 1 • • • > 4 5***5 lr / * 

One defines a suitable rational function w(t, z) with values in La [Aqo] as in |RVj . cf. 
|MV2llRSV| . The function is called canonical. The canonical function is symmetric 
with respect to the group £■ = S/ 1 x • • • x S* r of permutations of coordinates tj 
with the same upper index. One considers the system of equations 

(a n V' \- (a s ,ai) ^ (o^Ot) _» 

^•^ Z^ (i) + Z^ Z^ (i) _ ■» + Z^ (i) _ «) u ' 
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where i — 1, . . . , r and j = 1, . . . , li. The system is symmetric with respect to Sj. 
One shows that if t is a solution to (|6.1|1 . then w(t°,z) belongs to SingLA[Aoo] 
and w(t°, z) is an eigenvector of the Hamiltonians Ki(z), . . . , K n (z), see jRVj . 

This method of finding eigenvectors is called the Bethe Ansatz method. System 
(16.11) is called the Bethe Ansatz equation, the vector w(t° , z) is called a Bethe vector. 

The standard form of the Bethe Ansatz Conjecture says that if Ai, . . . , A n , Aoo 
are integral dominant, and z\, . . . , z n are generic, then the Bethe vectors form a 
basis in Sing La [Aoo]. In particular, the conjecture implies that the number of 
Si-orbits of solutions to l|6.1|l is not less than the dimension of Sing La [Aoo]- 

6.2. Specialized Bethe Ansatz equation. Fix a natural number r and a se- 
quence of non- negative integers m. Choose a sequence of non-negative integers I 
and a non-negative integer k. 

Consider the following system of algebraic equations for variables t — (t\^), 
i= j = l,...,r: 

(6.2) 

TUT + _ i ~ 2^ AT) TTy+Z^TT) T(3J~ U I* - V--i*iJ> 

l i l i 1 s, i ls s L i l s 

rrij \ - 2 x - 1 \ - 1 . . _ 

~Aj) ~ 1^ AS) _ Ai) + 2^ Jj) _ Aj+i) + 1^ Ai) _ (« = 

l i s, s=jti L i L s s l. t is s <-j ts 

77T " E (r) _ Ar) + E (r) _ Ar-1) = (* = 1, . . . , / r ), 

where j = 2, . . . ,r — 1. 

System Q6.2H is the Bethe Ansatz equation <|6.1[1 for the Gaudin model associated 
to the Lie algebra sL+i, points z\ = 0, — 1, and weights Ai, A2, such that the 
scalar products with simple roots are given by l|2.1|l . The Bethe Ansatz equation 
with two arbitrary points Z\ and z 2 is related to our choice of Z\ — and Z\ = 1 by 
a simple rescaling of variables t. 

System l|6.2f) is invariant with respect to the group X; = Ejj x • • • x S; r of 
permutations of variables with the same upper index. Thus the group acts on 
solutions of the system. This action is free and we do not distinguish between 
solutions in the same orbit. 

Assume that the weight 

Aoo = Ai + A 2 - ^2 l s a s 

s 

is dominant. Then it is known, see jMVlj . that the number of orbits of solutions of 
equation H6.2JI is finite and it is bounded from above by the dimension of the space 
of singular vectors of weight Aoo in the tensor product L\ ± ® La 2 . Using Lemma 
12.11 we obtain the following lemma. 

Lemma 6.1. If Aoo is dominant, then system l|6'.i| has no solutions unless m,l, k 
are consistent. In that case the system has at most one orbit of solutions. □ 

For t — {tl), we define a tuple of polynomials y = (yi, . . . , y r ), 

U 

s=l 



30 



E. MUKHIN AND A. VARCHENKO 



We are interested only in zeroes of yi and we consider each polynomial yi up to 
multiplication by a non-zero complex number. We say that y represents t. Note 
that degy 4 = k. 

We call a tuple of polynomials (jji ,y r ) generic if all yi have no multiple roots, 
no common roots with yi±\ and no common roots with Ti, where Ti are given by 
(|2.11|l . If y represents a solution t of the Bethe Ansatz equation then y is generic. 

6.3. Connection between Bethe Ansatz and spaces of polynomials. We fix 

a choice of consistent parameters m, I, k. Let U — U (m, I, k) and V = V(m, I, k) 
be the corresponding spaces of polynomials, see Lemmas 12.31 12.21 

Let fx,..., fr+i and gi,... ,g r +i be bases of V and U as in Lemmas 14.41 and |5"T1 
Then it is shown in |MV1| . cf. also Lemma [2.51 that polynomials W'(fi, ...,fi) 
and W'(gx, . . . , g r +x-i) are equal up to multiplication by a non-zero constant. Let 
yi be the monic constant multiple of W*(fx, . . . , fi). We have 

(6.3) Vl = cW^{h,...,h) (i=l,...,r), 

where c is a non-zero constant depending on i and fx, . . . , /j. 

For example in the case of lx = • ■ • = l r = 0, we have yi = ■ ■ ■ = y r = 1- 

Lemma 6.2. The Bethe Ansatz equation ()6'.i?|l for consistent m, I, k has a solution 
t if and only if tuple y — (yx, ■ ■ • ,Vr) is generic. Moreover, in such a case, y 
represents t. 

Proof. If the the Bethe Ansatz equation l|6.2|l . associated with consistent parame- 
ters m,£, k, has a solution t, then the solution t generates "a population of solu- 
tions", see Section 3.4 in |M V1| . That population of solutions has its "fundamental 
space" Vt, which is an r + 1-dimensional vector space of polynomials with exactly 
two singular points at and 1 and exponents given by (|2.6(l . see Section 5.3 in 
|MV1| . Thus Vt — V(m,l,k) by Lemma [2.21 Having the fundamental space Vt, 
define the tuple y by formula l|6.3l) . Then y represents the solution t, see Theorem 
5.12 in |MV1| . Since t is a solution of the Bethe Ansatz equation, the tuple y is 
generic in the sense of Section 16.21 This proves the lemma. □ 

6.4. Properties of polynomials yi, . . . ,y r . Let m, I, k be consistent parameters. 
We have the corresponding monic polynomials yi (m,l, k), see (|6.3|l and Lemma 

m 

Proposition 6.3. For all consistent parameters m, I, k we have 

yx(m,l,k)(l)^0, yi (m,l,k)(0)^0 (i = l,...,r). 

Proof. The proposition holds if and only if F r+ x- s (0) D F s (oo) — and F r (l) n 
^1(00) = in V(m,l,k). The proposition is obviously true if (I, k) — (0,0). Now 
the proposition follows for all consistent parameters m, I, k because there exists a 
canonical map V(m, 0, 0) — > V(m, I, k). □ 

Note that Proposition l6 . 31 shows that polynomials yi and Ti are relatively prime. 

Lemma 6.4. The polynomials yi(m, l,k), i = 1, . . . , r, are polynomials of degree 
li whose coefficients are rational functions of m,k with coefficients in Q. 

Proof. From Theorem 13.31 it is clear that we have a basis in V(m,l,k) which is 
compatible with the flag at infinity and whose coefficients are rational functions of 
m, k. Clearly, we still have rational coefficients after the operation of taking the 
divided Wronskian. □ 
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Lemma lfi.4l allows us to consider polynomials yi(m,l, k) for all generic complex 
parameters m, k. 

Next two lemmas describe relations between polynomials corresponding to dif- 
ferent r. 

Let fh = (m 2 , . . . , m r ), I = (l 2 , ■ ■ ■ , l r )- 
Lemma 6.5. We have 

lim yi(m, I, k) = (x — 1) , 

mi — >oo 

lim yi(m,l,k) = yi-xtfhjjx), 

mi — >oo 

where i = 2, . . . , r. 

Proof. Recall the basis of V(m,l,k) described in Lemma [4.41 From l|4.2(l we see 
that Vi have finite limits when m\ tends to infinity. Denote these limits by Vi. From 
Proposition ^. 101 we conclude that vq = (x — 1) 1 . Consider the space 

{vi, x m2+1 v 2 , . . . , a£*=> m *+ 4 - V}- 

This space has the same exponents as V(fh, Indeed, it is clear about expo- 

nents at infinity. The exponents in the finite points can only increase in the limit. 
It does not happen because of relation (|2.5|l . 

Therefore this space coincides with V(fh,l, li). The lemma follows. □ 

Let rh = (mi, . . . , m r , n), where n is any complex number, I — (li, . . . , Z r , 0). 
Lemma 6.6. We have 

yi(rh,l,k) = yi(m,l,h), 
y r+1 (m,l,k) = 1, 

where i = 1, . . . , r. 

Proof. The lemma follows from Lemma 13.131 □ 
Finally we describe properties of zeros of y{ in a certain asymptotic zone. 

Lemma 6.7. Fix a partition I and a negative real number k. Let m be a sequence 
of real numbers such that 

mi << ni2 << ■ ■ ■ << m r << 0. 

Then the corresponding tuple y is generic. Moreover all yi have real distinct roots 
t^ such that < t [ p < 1 and t^ > t^ ] if j x < j 2 . 

Proof. If all rrii < then all roots of yi(m,l,k) are distinct, real and located in 
(0, 1). This is a general fact about roots of multiple orthogonal polynomials, see 

ch]. 

Let mi << 0, and let all other mi to be real. Then the roots of y 2 tend to 
the roots of another multiple orthogonal polynomial by Lemma lr!"51 Therefore the 
roots of y 2 are also real and distinct. Moreover, they are smaller then the roots of 
2/i, which all are close to one. 

Let mi << m 2 << 0, and let all other m, to be real. Then the roots of j/3 also 
tend to the roots of a multiple orthogonal polynomial by Lemma [6.51 Therefore 
the roots of j/3 are also real and distinct. Moreover, they are smaller then the roots 
of y 2 , which are close to one in this zone. 

The lemma follows. □ 
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6.5. Bethe vectors. Let h, . . . , l r , k be non-negative integers such that k > h > 
■ ■ ■ > l r . Let rrii be real numbers. 

We list properties of the canonical function used for the construction of the Bethe 
vectors w(t) = w(t,0, 1) e £a[Aoo], see |RSV| . 

(i) We have i A [Aoo] = ffi w + P2 =A C o i A 1 M ®£a 2 [M2]- For evei T sucn Mi>M2, 
there exist a rational function w^ 14L2 {t) £ ® L\ 2 [/i 2 ] such that 

«'(*) = E w+Ma =A««'w,/* a (*)- 

(ii) Let Fx, • • • , F r 6 tt_ be the standard generators. Let f7(n_) be the complex 
universal enveloping algebra of n_. Let £/ R (n_) C U{xi-) be its real part 
(consisting of polynomials in Fx, . . . ,F r with real coefficients). For s = 1,2, 
let v s 6 L\ s be a highest weight vector. Let L R be the real part of La s , 
i.e. L R is the image of v s under the action of C/ R (n_). For all /ii,/i2, we 

have tu^j i/J2 (i) £ [/ix] <8> £ R 2 [^2] if ah coordinates tj of the vector t are 
real. 

(iii) Let = be the product of symmetric groups. Let v be a function of 
£ . We define the action of n 6 on ^ by permuting the ij*' 's with the same 
upper index. We define the symmetrizer operator Sym v = E 7reS; tt • v. 
We have the following formula 

v(t) v x ® Ft- h [F 2 ,F 1 } 1 ^ . . . [F r , [F r _x, [. . . , [F 2) Fx] . . . 



v{t) 



where 

1 



{h-h)\ ... (Jr-l-ir)! M 
^1 — ^2 -1 ^2 — 

n 7i(T) tt n 



n 



1 (\ J-J i=X l I i+ ; 1 _; 2 J-J I 1 ! l i+l 1 -l 2 ) 
1 



^^1 -lVi (2) )(t {3) -t (2) 1 ftM - 1^- 1 ) } J ' 

Theorem 6.8. Fix a partition I. Then for generic complex values of k, mi, . . . , m r 

we have: 

• system (|6'.2[) has a unique orbit of solutions, 

• the multiplicity of this solution is equal to one, 

• the corresponding Bethe vector is non-zero. 

Proof. Consider the tuple y(m, I, k) — (yi(m, I, k), . . . , y r {fn, I, k)) associated with 
parameters m,l, k. By Lemma [6.71 the tuple is generic for generic complex values 
of the parameters m, k. Hence for almost all consistent values of parameters m, I, k 
with fixed I, the tuple y{m, I, k) is generic. The exceptions lie in a proper algebraic 
subset in the space of all m, k. 

Hence by Lemma 6.1 for almost all consistent m,l,k with fixed I the Bethe 
Ansatz equation has exactly one orbit of solutions. Equisingularity theorems in 
[Vlj imply that the Bethe Ansatz equation has exactly one orbit for almost all 
complex values of parameters m, k with fixed I. 

Let us prove that the solutions of that only orbit have multiplicity one. Indeed 
if m,l,k are as in Lemma 16.71 then using the properties of the weight function 
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described above, we obtain that the corresponding Bethe vector w(t) is non-zero 
and belongs to the real part of La- Then the multiplicity of solutions of the only 
orbit is equal to one by Corollary 7.2 in |MV2| . Hence the multiplicity of solutions 
of the only orbit is equal to one for almost all values of m, k. On the other hand if 
the multiplicity is one, then the corresponding Bethe vector is non-zero, see MV2 . 

□ 

6.6. A counterexample to the Bethe Ansatz Conjecture. According to The- 
orem f° r generic values of parameters, there is a unique solution of the Bethe 
Ansatz equation (|6.2|l and a unique non-zero Bethe vector. In particular, for generic 
values of m, k the Bethe Ansatz Conjecture holds. In this section we give an ex- 
ample of parameters when this is not true. 

We consider the case of r = 2, l\ = 2, I2 = 1. Then polynomials y\ and y 2 are 
computed to be 

_ m 2 (mi+m 2 + l) 

y2 (m 2 + 2)(k + mi + m 2 ) 

and 

yi = (x 1} 2 H (fc-l)(2fc + m 2 + 2m 2 ) fc(fc - 1) 

(mi + k — l)(mi + m 2 + k) (m 2 + mi + k)(mi + k — 1) 

Lemma 6.9. Suppose (2mi + m 2 ) 2 + fc(4mi — m|) = 0. Then we have y\ = y\. 
In particular the pair (yi,y 2 ) is not generic. □ 

Let 

r = 2, k = 49, mi = 2, m 2 = 3, Zi = 2, l 2 = 1. 
Define weights Ai, A 2 by the conditions (Ai,a») = m^, (A 2 ,ai) = kSu. Let 
= Ai + A 2 — liai — Z 2 a 2 . Let L\ be the irreducible sl^ module of highest weight 

A. 

Proposition 6.10. The multiplicity of L\ x in L\ t ®L^ 2 is 1. The corresponding 
Bethe equation l|6'.ij|) has no solutions. 

Proof. The parameters m,l,k are consistent, therefore the multiplicity of La^ in 
Lai &> La 2 is one by Lemma l2~Tl The pair yi(m, I, k), y 2 (m, Z, fc) is not generic by 
Lemma fo. 91 therefore system (|6.2|l has no solutions. □ 

Proposition l6.10l provides a counterexample for the standard version of the Bethe 
Ansatz Conjecture for the Gaudin model. 

Remark. Here is another counterexample to the standard version of the Bethe 
Ansatz conjecture. 

Consider the tensor product of two adjoint sl^ representations. Then the mul- 
tiplicity of the trivial module in this tensor product is 1. One can explicitly check 
that the corresponding 3-dimcnsional vector space V is 

\Z = span{(2a;-l) 2 , x 4 }. 

The space V is the unique 3-dimensional space of polynomials spanned by polyno- 
mials of degrees 2, 3 and 4 which has two finite singular points 0, 1 such that the 
exponents at both singular points are 0, 2,4. 

Then the corresponding pair (yi,y 2 ) = ((x — l/2) 2 ,(x — 1/2) 2 ) is not generic 
and by the same argument as in Proposition 16 . 10l the corresponding Bethe Ansatz 
equation has no solutions. 
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7. Appendix 

7.1. A scalar product. We compute the scalar product of yi and y r with respect 
to the measure defined by integration with the Jacobi type weight. 
Let 

tJ m ,k = (x-l)- k - 1 x- m -\ 

Set 

r 

\m\ = to, + r — 1. 

i=i 

Then we have wi m |,fc = w, where u> is given by ^4.1[l . 

Assume that m and fc are negative real numbers. Define a scalar product of 
functions f(x) and g{x) by the formula 

(f(x),g(x)) m .k = / f{x)g{x)u m>k {x)dx. 
Jo 

The Jacobi-Pineiro polynomial P(m, I, k) is the unique monic polynomial of degree 
li which has the property 

(P,O| m |, fe = 

for all s of the form 

r 

s= ^2 rn,j + t, i e {0, 1, . ..,r — 1}, t — i e {0, 1, .. .,l r -i - l r -i+x — 1}, 

j—r—i+l 

where k is a negative number and rrij are such that Y] -_ mj + i — 1 < for 
i= 1, ...,r. 

Operators Di and are adjoint in the following sense: 
Lemma 7.1. We have 

(Di(m,l,k)f,g)\ m \ +1>k+1 = -{f,D^(m,l + l h k + l).g)| m | + i, fc . 
Proof. The lemma is checked by integration by parts. □ 

Using Lemma 17. II and Propositions 14.101 and 15.41 the following scalar product is 
easy to compute explicitly: 

Lemma 7.2. We have 

r-l 

(yi(m,l,k),y r (m,l,k))\ m \ +ltk = (1, l)\ m \+i,k-h x II 

i=0 

TT Lt=r-j m * + 1 - 3 TT Et=i+l TO *+ r + 1 -' 

11 , . i . / , / 11 



. 3- 

□ 
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7.2. Afflne Weyl group. In this section we note that operators Di together with 
the Weyl group generate an extension of the affine Weyl group. 

We call parameters m, I, k admissible if the polynomials Vi(rn, I, k), i = 0, . . . , r, 
are well defined. For admissible parameters m, I, k, we define the corresponding 
space V(m,l,k) as follows. The space V(m, I, k) is the space of functions with 
the explicit linear basis {vo(m, I, k), vx(m, I, k)x mi+1 , . . . , v r (m, I, k)x^i =1 m i+ r } ) 
where Vi(m, I, k) are given by 14.2(1 . 

Let V be a disjoint direct sum of spaces V(m, I, k), 

V = ® m ,i,kV(m,l,k) 

where the sum is taken over all sequences m of r complex numbers, all sequences 
I of non-negative integers and all non-negative integers k such that m, I, k are 
admissible and k > 1% > ■ ■ ■ > l r . 

The Weyl group S r +i acts on V by acting on parameters. Namely, we define the 
action of the Weyl group on V by the following rule: 

Si : V(m,l,k) —>V(si-m,(si) k l,k), 
Vj (m, I, k)x^l=r m *+3 i > v Si{j) ( Si • rn. ( Si ) k l, k)x^*=" 

where i = 1, . . . , r and j = 0, . . . , r. 

Introduce linear operators Di, i = 0, . . . , r, acting on V. We set 

Di : V(m,l,k) -> V(m,l + U,k+ 1) 

so that the restriction of operator Di to V(m, I, k) equals to Di(m, I, k). 
We have the following commutation relations. 

Lemma 7.3. We have 

(7.1) DiDj - DjDi, SiDj - D s . u)Si . 

□ 

Consider the semigroup G generated by Si, i — 1, . . . , r, and Di, i = 0, . . . ,r, 
subject to the Weyl group relations s| = id, SiSi+iSi — Si+iSiSi+%, SiSj — SjSi 
(I* — j\ > 1) and relations Q7.1J1 . 

Note that the element Do . . . D r is in the center of G. The factor of G by the 
sub-semigroup generated by this element is isomorphic to the affine Weyl group of 
sl r+ i, cf. |K) . Proposition 6.5. 
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